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ABSTRACT 

Lower bounds are derived on the amplitude B of intergalactic magnetic fields (IGMFs) in the 
region between our Galaxy and the blazar Mrk 421, from constraints on the delayed GeV flux 
of pair echos that are emitted by secondary e~e^ produced in 77 interactions between primary 
TeV gamma-rays and the cosmic infrared background. The distribution of galaxies mapped by 
the Sloan Digital Sky Survey shows that this region is dominated by a large intergalactic void. 
We utilize data from long-term, simultaneous GeV- TeV observations by the Fermi Large Area 
Telescope and the ARGO-YBJ experiment extending over 600 days. For an assumed value of B, 
we evaluate the daily GeV flux of the pair echo expected from the TeV data, select the dates where 
this exceeds the Fermi 2-a sensitivity, compute the probability that this flux is excluded by the 
Fermi data for each date, and then combine the probabilities using the inverse normal method. 
Consequently, we exclude B < lO"^*^ *^ G for a field coherence length of 1 kpc at ^ 4-ct level. This 
is much more significant than the 2-a bounds we obtained previously from observations of Mrk 
501 involving Cherenkov telescopes, by virtue of more extensive data and improved statistical 
analysis. Compared with most other studies of IGMF bounds, the evidence we present here for a 
non-zero IGMF is more robust as it does not rely on unproven assumptions on the primary TeV 
emission during unobserved periods. 

Subject headings: magnetic fields — gamma rays: observations — galaxies: active — gamma rays: 
theory — BL Lacertae objects: individual (Mrk 421) — radiation mechanisms: nonthermal 



1. Introduction 

Intergalactic magnetic fields (IGMF), espe- 
cially those inside intergalactic void regions, 
have attracted much attention because they 
could be a remnant of phenomena that gener- 
ated prim ordial magnetic fields in the early Uni- 
verse (e.g.lGnedin et aL 200dlLanger et al. 2005t 



Takahashi et al. 20051: llchiki et al. 20061 ). While 
primordial magnetic fields can be amplified later 
within galaxies and galaxy clusters by dynamo 
processes, they may remain unaffected by sub- 
sequent astrophysical effects deep inside voids. 
Thus, IGMFs are expected to be a window onto 
the early Universe. For comprehensive reviews on 
primordial and intergalactic magnetic fields, see 
Widrow (2002), Widrow et al. (2012) and Ryu et 
al. (2012). 

However, the predicted amplitudes of primor- 



dial magnetic fields are generally very small, 
B = 10^25 _ 20^15 Gauss, and difficuh to probe 
with methods such as Faraday rotation of dis- 
tant radio sources or the anisotropy of cosmic 
microwave background. In this context, a novel 
method utilizing delayed secondary emission from 
high-energy gamma-ray sources was proposed by 
Plaga (19 95) and subseque n tly developed by many 

authors JPai et al. 20021: iRazzaaue et al. 2004t 

Murase et al. 20071 ' Murase et al. 2008Hlchiki et al. 2008t 



Takahashi et al. 



Elviv et al. 2009HNeronov fc Semikoz 2009t 



Takahashi et al. 20111 ). This emission that we re- 



fer to as "pair echo" is expected to occur typically 
at GeV energies, for which the Fermi Large Area 
Telescope (LAT) provides important constraints. 
Since the echo flux is predicted to be larger for 
smaller B, a GeV upper limit translates into a 
lower bound on B. 



In our previous study ([Takahashi et al. 20121) . 
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Fig. 1. — Distribution of galaxies between Mrk 
421 and our Galaxy mapped by SDSS. The pluses 
and crosses represent galaxies located within a 1- 
degree slice and a 3-degree slice (but outside the 1- 
degree slice) from the plane with declination 38.2, 
respectively. The width of 1-degree slice is about 
1.7 Mpc at 100 Mpc from our Galaxy. 



we focused on a specific TeV flare of Mrk 501 
observed in 2009 by the VERITAS (Very En- 
ergetic Radiation Imaging Telescope Array Sys- 
tem) and MAGIC (Major Atmospheric Gamma- 
ray Imaging Cherenkov) telescopes. Calculating 
the lightcurve of the pair echo expected from the 
flare and the quiescent emission before and af- 
ter the flare and comparing with simultaneous 
Fermi observations, we obtained a lower bound 
on the IGMF of B > lO^^o Gauss (90% C.L.) for 
a field coherence length of 1 kpc. Being based 
firmly on the observational data and nearly free 
of assumptions concerning the primary TeV flux 
in unobserved periods or spectral bands, our re- 
sult is more robust com p ared to previous studies 
(iNeronov fc Vovk 2010l: lAndo fc Kusenko 20l"ot 



Dermer et al. 2011tlNeronov et al. 2011tlTavlor et al 



IGMFs. Mrk 421 has been monitored continu- 
ously at TeV energies by the ARGO-YBJ exper- 
iment over the period from 2007 November to 
2010 February ([Bartoli et al. 201 II ). during which 
many flares were observed so that more statis- 
tically significant bounds on IGMFs can be ex- 
pected. Note that compared to Cherenkov tele- 
scopes, such air shower detectors have a much 
higher duty cycle and allow uninterrupted long- 
term observations, albeit at lower sensitivity. 

2. TeV and GeV Emission from Mrk 421 

For the calculation of pair echo, we need the 
spectrum and lightcurve of TeV emission. In Bar- 
toli et al. (2011), the daily fluxes above 0.3 TeV 
are given for about 1000 days. Among them some 
negative daily fluxes are reported. They would be 
due to the systematic error and we simply set them 
zero. Although the daily spectra are not given, 
the average spectra were derived from four flux- 
level data groups according to the X-ray counting 
rate. In this letter, because the X-ray data is not 
available, we devide the data into three groups ac- 
cording to the TeV daily count. This is justifled 
because the TeV count was shown to be tightly 
correlated with the X-ray count rate. Our three 
groups, "high" (TeV daily count > 80), "medium" 
(80 > TeV daily count > 40) and "low" (TeV daily 
count < 40), correspond to the X-ray flux levels 4, 
3 and 2 (level 1 has almost the same spectral in- 
dex with that of level 2 so low state includes level 
1 as well), respectively. We adopt the spectral in- 
dices given by Bartoli et al. (2011) and scale the 
normalization with the flux day by day according 
to the corresponding group. Minimum and max- 
imum cutoffs are also imposed at 0.1 TeV and 5 
TeV, respectively, the latter corresponding to the 
highest energy photons detected by ARGO-YBJ. 



Tavecchio et al. 2011l:lTavecchio et al. 20101: iDolag et al. 3&lHt g. El absorbed and unabsorbed spectra of the 

201lt iree groups are shown. Here the absorbed spectra 



Arlen et al. 2012f ). 

Here we focus on the TeV blazar Mrk 421 lo- 
cated at z = 0.031. As seen in Fig. [U maps of the 
local galaxy distribution from the Sloan Digital 
Sky Survey reveal that a huge void lies between 
our galaxy and the superc l uster containing Mr k 
421 (jAbazajian et al. 20091 : iBlanton et al. 20051 ). 
This is also seen to be the case for Mrk 501. 
Thus, Mrk 421 is a desirable target for probing 



are calculated using the CIB model of Franceschini 
et al. (2008), which was also adopted in Bartoli et 
al. (2011). It turns out that most of constraints on 
IGMF come from the emission during high state 
and low state does not affect the result at all. 

On the other hand, concerning GeV gamma- 
rays, Fermi LAT has performed uninterrupted 
monitoring of Mrk 421 in the survey mode from 
MJD 54683. The data were acquired via the Fermi 
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Fig. 2. — The unabsorbed (solid) and absorbed 
(dashed) spectra of Mrk 421 for the high, medium 
and low states (from top to bottom). The values 
in the parentheses indicate the daily TeV count. 

Science Support Center (FSSC). We use the stan- 
dard analysis tools supplied by FSSC in three en- 
ergy bands, 100 MeV-1 GeV, 1-10 GeV, and >10 
GeV, so as to keep reasonable photon statistics 
in time intervals as short as one day while retain- 
ing some energy resolution. Since the statistics is 
still small, we adopted the aperture photometry 
method, i.e., events falling within one degree from 
the source were counted. Note that background 
events above 1 GeV in one-day bins at the high 
Galactic latitude of Mrk 421 are essentially negli- 
gible. Flux probability distribution functions were 
calculated assuming Poisson statistics when the 
data show no significant gamma-ray signal. Below 
we use data during MJD 54683 - 55255, when both 
TeV and GeV observations were performed. Also 
we focus on an energy range of 1 — 10 GeV for the 
pair echo, where Fermi LAT is the most sensitive 
to the pair echo. 

3. Pair Echo 

Here let us summarize briefly the basic physics 
of pair echo emission (see e.g. Ichiki et al. 
2008 and Takahashi et al. 2012 for more de- 
tails). Primary gamma-rays with energy E-y > 
1 TeV emitted from an extragalactic source 
have mean free path A-,^ — l/(0.26crT7iiR) = 
190 Mpc (riiR/0.01 cm^'^)^^ for 77 pair pro- 
duction interactions with photons of the cos- 
mic infrared background (GIB), where ctt is the 
Thomson cross section and njR is the number 



density of CIB photons most relevant for the 
interactions. The produced pairs with energy 
Ef. « £^7/2 give rise to the pair echo emission 
by IC upscattering of ambient cosmic microwave 
background (GMB) photons to average energy 
(E^echo) = 2.7rcMB7e = 2.5 GcV {Ey/2 TeV)2, 
where 7e = Ee/meC^ is the Lorentz factor of the 
pairs and Tcmb = 2.7 K is the GMB temper- 
ature. Thus, primary gamma-rays in the range 
Ej ~ 1 — 5 TeV induce echos with typical en- 
ergies Eccho — 1 — 10 GeV. The pairs upscat- 
ter GMB photons successively until they lose 
most of their energy after propagating an IC 
cooling length Aic,cooi = ^rnl/{AEeCrTUcMB) = 
350 kpc {Ee/l TeV)-i, where Ucmb is the GMB 
energy density. The length scales for A^^ and 
Aic,cooi imply that the secondary pairs typically 
arise in locations far removed from the source 
on scales of intergalactic voids, whereas the pairs 
propagate only for short distances within such re- 
gions while generating the echo emission. In our 
case with Mrk 421, as we saw in Fig. [U there is a 
huge void between our galaxy and Mrk 421 and, 
noting that the typical scattering angle of gamma- 
rays and charged particles is ^ l/7e ^ 10~^, their 
propagation region is mostly in the void space. 

The pair echo emission arrives at the observer 
with a time delay relative to the primary emis- 
sion, caused by the effects of angular spreading 
in pair production and IG interactions, as well 
as by deflections of the pairs in intervening mag- 
netic fields. The typical delay time due to angu- 
lar spreading is Atang = (A-y^, -I- Aic,cooi)/27e « 
3 X 10^ sec (-Eocho/1 GeV)-H'^iR/0-01 cm-^)-! 
(Ichiki et al. 2008), while that due to magnetic 
deflections is Ate — (A^^. -|- Aic,cooi)(^'b)/2, where 
(6»|)V2 = max[Aicxooi/?'L, (Aic,cooircoh/6)^/V^L] 
is the variance of the magnetic deflection angle in 
terms of the Larmor radius tl and field coherence 
length Tcoh- If Tcoh < Aiccooi, 

AtB « 2 X lO'^ sec (£;ccho/l GeV)-^/^{B/10-^^ G)^ 
x(reoh/l kpc)(niR/0.01 cm-3)-i, ^-^^ 

where B is the field amplitude. Hereafter we fidu- 
cially take rcoh = 1 kpc (see e.g. Langer et al. 
2005), but the results are scalable to other values 
of Tcoh since it enters only through the combina- 
tion B^rcoh and only when rcoh ^ Aic,cooi- The to- 
tal delay time is approximately At — A^ang + A^b , 
and the magnetic field properties are refiected in 
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the delay as long as A^ang ^ Ate- 

The spectra and light curves of the pair echo 
can be evaluated as follows. For a primary fluence 
dN^/dEj, the associated time-integrated flux of 
secondary pairs is 



-(7ej = 4me— -(i/^, = 2TOe7e) 



d-je 



(2) 



where T~f^{E^) is the 77 optical depth in the CIB. 
The time-dependent echo spectrum is 



echo 



dtdE^ 



djc 



dN, d^Nic 
dje dtdEj ' 



(3) 



where d^Nic/dtdE-y is the IC power from a sin- 
gle electron or positron, and dN^jd^e is the to- 
tal time- integrated flux of pairs responsible for the 
echo emission observed at time related nontriv- 
ially to dNf, Q/d^e in Eq. ©. In Takahashi et 
al. (2012), we developed a formalism to calculate 
this factor accounting for the finite probability of 
pair production near the observer, extending the 
earlier work by Dai et al. (2002). 

Note that the pair echo fluence is determined 
by the total amount of absorbed primary gamma 
rays and thus independent of the IGMF, in con- 
trast to the pair-echo flux which is roughly given 
by the fluence divided by At. Weaker IGMFs 
generally give higher echo fluxes, as long as the 
time delay does not become dominated by angu- 
lar spreading and the echo flux remains sensitive 
to B. For Tcoh = 1 kpc, A^b approaches Atang 
if i? 10^^*^ G. This behavior is clearly seen in 
Fig. [3] which shows the lightcurves of pair echo 
with 1 — 10 GeV and TeV daily counts. Pair-echo 
flux rises on the occasion of TeV emission and de- 
cays in the absence of it. For smaller IGMF the 
response to TeV emission is quicker and the peak 
flux is larger. Finally, here, we do not consider the 
possibility that the pairs lose energy by heating 
the intergalactic gas t hrough plasma instabilities 
(jBroderick et al. 20121 ). as the efficiency of this 
process has been disputed ([Miniati k. Elviv 20121 ) 
and is highly uncertain. 

4. Statistical Analysis 



In this section, we compare the calculated pair- 
echo flux with the Fermi-LKY data and derive 
constraints on IGMF. In contrast to the previous 
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Fig. 3. — Expected lightcurves of the pair echo 
in the 1—10 GeV band (solid and dashed lines 
for B = lO"^"'^, 10"^^'^ Gauss, respectively) com- 
pared with the observed TeV lightcurve (points, 
arbitrary units). 



paper ( Takahashi et al. 2012[ ). we have vast num- 
bers of independent flux bins (each bin represents 
the daily flux) so we need a more sophisticated 
way to derive constraints. First, we compute the 
probability Pi that a specific value of IGMF am- 
plitude is excluded by the i-th flux bin, using the 
probability distribution function of the true flux 
obtained from Fermi-LAT observation. Then, we 
combine the probabilities to derive the total prob- 
ability Ptot using meta-analysis. 

Here, it should be noted that it is not appro- 
priate to combine probabilities from all bins. The 
reason is as follows. If the TeV emission is very 
weak for the i-th bin, the expected pair-echo flux 
on that bin is very low irrespective of the value 
of IGMF. Consequently, given an upper bound on 
GeV flux by Fermi, the probability Pi would be- 
come very small. If we combine such P^s, the to- 
tal probability Ptot would also become very small, 
that is, any values of IGMF would not be excluded 
even if some P^s are very large at TeV flares. In 
fact, it is not fair to combine P^s of time bins 
where the pair-echo fluxes are expected to be much 
smaller than the Fermi sensitivity because it is 
natural for the pair echo not to be observed and 
Fermi data would not have any information on 
IGMF. Thus, we need to choose data bins where 
the pair echo should have been observed. Note 
that the data bins which should be chosen depends 
on the value of IGMF because the pair-echo flux 
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Fig. 4. — Fermi-hAT 2-a daily sensitivity 
(dotted), lightcurves of pair echo with B — 
10-20-5 Gauss (soHd) and B = lO^^o Gauss 
(dashed) and Fermi 50% upper Hmits (points). Ah 
are for an energy range of 1 — 10 GeV. 



depends on it. In general, as we explained above, 
larger IGMF results in a lower flux so that the pair 
echo is harder to be observed and the number of 
chosen data bins would be smaller. 

In this letter, we set a threshold to choose data 
bins such that the expected pair-echo flux exceeds 
the 2-cr sensitivity of Fermi-L AT. In Fig. gl the 2- 
(T sensitivity and the pair-echo lightcurves for B = 
10-20-5 Gauss and B = IQ-^o Gauss are plotted 
for an energy range of 1 — 10 GeV. As can be seen, 
during this period (which is only a small part of 
the whole period), 4 and 3 bins exceed the Fermi- 
LAT sensitivity in cases with B = lO-^O-S Gauss 
and B = lO-^o Gauss, respectively. Looking back 
at Fig. [31 large TeV flares can be seen at these 
selected bins. 

In Fig. m the 50% upper limits on the daily 
flux are also plotted. Looking at the first flare 
at MJD 55147, the Fermi-LAT upper limit is 
much smaller than the expected pair-echo flux for 
B = 10-20-5 Gauss so that the probability that 
this value of IGMF is excluded is very large. On 
the other hand, for B — lO^^o Gauss, the pair- 
echo flux does not exceed the sensitivity and this 
bin is not counted to compute Ptot, even though 
the pair-echo flux is larger than the upper limit. 
As to the second (MJD 55152) and third (MJD 
55166) flares, the pair-echo fluxes are larger than 
the upper limits both for B = 10^20-5 Gauss and 
10-20 Gauss. At the fourth flare (MJD 55182), 



the pair-echo flux for B — 10-20-5 Gauss is almost 
the same as the upper limit. Because this is the 
50% upper limit, this bin does not favor or ex- 
clude this value of IGMF. On the other hand, in 
the case with B = IO-20 Gauss, the pair-echo flux 
is smaller than the upper limit and this value of 
IGMF is allowed. 

In this way, for a specific value of IGMF, appro- 
priate bins are chosen and the expected pair-echo 
flux is compared with the Fermi-LAT data. Here, 
we consider the probability distribution function 
of the true flux and calculate the probability. Pi , 
that the true flux is below the pair-echo flux for 
the i-th bin. Then, to combine P^s, we use the 
inverse normal method, which is a kind of meta- 
analysis as we explain below. First, derive the Z 
value of the normal distribution for the i-th bin, 
Zi, which corresponds to the percentile (point) of 
the one-sided P value Pi. It should be noted that 
Zi is negative if Pi < 0.5. Second, compute the 
total Z value Ztot by the following equation. 



Ztnt — 



1 ^ 



(4) 



where iV is the number of the selected bins. Fi- 
nally, derive the one-sided P value Ptot of the nor- 
mal distribution which corresponds to the above 
Ztot- Thus, we can interpret this Ptot such that 
this value of IGMF is excluded at the confidence 
level Ptot- 

Fig. E] shows the Z value Ztot as a function of 
the IGMF amplitude. Because the time delay At 
is dominated by the angular spreading time for 
B < 10-20-5 Gauss, the result does not change 
for that range and such a weak IGMF including 
zero IGMF is excluded by about 4-cr significance. 
The significance decreases for larger IGMFs and a 
IGMF larger than 10-^^'^ Gauss is not constrained 
at all. This is because the pair-echo flux is strongly 
dependent on the IGMF value and there are no 
bins where the pair-echo flux exceeds the threshold 
for B > 10-19-5 Gauss. 

Finally, it is noticed that we assumed there is no 
emission in GeV range other than the pair echo. In 
fact, there must be the primary GeV emission and 
possively some other emission. If we could know 
their existence and the amount somehow, the con- 
straints on the pair echo would become tighter, 
given the same upper limit on the GeV emission 
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Fig. 5.— The Z value Ztot that the value of IGMF 
is excluded by the Fermi data. 

by f ermi-LAT. In this letter, we consider only pair 
echo for GeV emission, which leads to conservative 
constraints on the amount of pair echo and then 
IGMF. 

5. Discussion and Summary 

Using data from simultaneous GeV-TeV obser- 
vations of Mrk 421 by fermi-LAT and ARGO- 
YBJ, we constrained the GeV flux of pair echos 
and derived lower bounds on the IGMF strength 
inside a large void region located between our 
Galaxy and Mrk 421. This was done by: 1) calcu- 
lating the daily pair-echo flux from the TeV data 
over 600 days, 2) selecting the dates where the ex- 
pected pair-echo flux exceeds the Fermi-LKT 2-cr 
sensitivity, 3) computing the probability that an 
assumed value of the IGMF is excluded by the 
i^ermi-LAT data for each date, and 4) combining 
these probabilities to derive the total probability 
using the inverse normal method. Consequently, 
an IGMF weaker than 10"^'^-^ Gauss for a field co- 
herence length of 1 kpc is excluded by about A-a. 
The result can be roughly scaled for other values of 
r-coh as B > 10-22 max[(rcoh/350 kpc)-!/^ 1] G. 

As with our previou s analysis using Mrk 
501 (|Takahashi et al. 20121 ) . no assumptions are 



made here concerning the TeV emission during 
unobserved periods. The obtained constraints are 
thus more robust compared to other studies, in 
particular those based on limits to the spatially- 
extended halo emission from secondary pairs that 
inevitably involves very long time delays, often 



longer than the typical lifetimes of blazars. Al- 
though the value of the lower limit obtained here 
is similar to our previous work, the statistical sig- 
nificance has improved remarkably, from less than 
2-cr to about 4-(T, thanks to the much larger data 
and improved statistical analysis. 

Here we used the Fermi-LAT data only as upper 
limits to the GeV fluxes. Because the pair-echo 
flux is strongly dependent on the TeV flux, we 
can potentially obtain tighter constraints on the 
IGMF by searching for statistical correlations be- 
tween the Fermi-hKI data and ARGO-YBJ data. 
This will be presented elsewhere in near future. 
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